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An  atomistic  simulation  technique  is  performed  to  investigate  the  molecular  structure  and  transport 
dynamics  inside  a  hydrated  Nation  membrane  and  a  hydrated  sulfonated  poly( ether  ether  ketone 
ketone)  (SPEEKK)  membrane.  The  simulation  system  consists  of  the  representative  fragments  of  the 
polymer  electrolytes,  hydronium  ions  and  solvent  molecules,  such  as  water  plus  methanol  molecules. 
Simulation  results  show  that  the  hydrated  SPEEKK  has  less  phase  separation  among  hydrophobic  and 
hydrophilic  regions  in  comparison  with  the  Nation.  Those  water  channels  formed  in  the  SPEEKK  are 
much  narrower  compared  to  those  in  the  Nation.  These  characteristics  lead  to  a  lower  mobility  of 
hydronium  ions  and  water  molecules  and  hence  relatively  lower  diffusion  coefficient  of  methanol  in 
the  SPEEKK.  It  results  in  the  reduction  of  the  methanol  permeation  problem  in  direct  methanol  fuel 
cells. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Polymer  electrolyte  based  fuel  cells,  including  proton  exchange 
membrane  fuel  cells  (PEMFCs)  and  direct  methanol  fuel  cells 
(DMFCs),  are  considered  as  the  most  promising  power  sources 
for  portable  electronic  devices  and  future  clean  vehicles  [1].  Cur¬ 
rently,  Nation  membranes  manufactured  by  Dupont  are  the  most 
popular  electrolytes  used  in  PEMFCs  and  DMFCs.  This  is  because 
of  the  acceptable  proton  conductivity  as  well  as  high  chemi¬ 
cal  and  mechanical  stabilities.  However,  high  cost  and  serious 
methanol  permeation  have  limited  Nation  and  other  similar  per- 
fluorinated  membranes  for  commercial  application  in  DMFCs. 
Over  the  past  decade,  there  were  many  studies  to  develop  non- 
perfluorinated  polymers  as  the  alternative  membranes,  such  as 
sulfonated  poly(aryl  ether  ketone)s,  sulfonated  polysulfones,  sul¬ 
fonated  polybenzimidazoles,  sulfonated  polyimides  and  some 
others  [2].  Sulfonated  poly( ether  ether  ketone  ketone)  (SPEEKK)  is 
a  member  of  the  sulfonated  poly(aryl  ether  ketone)  family  and  is 
regarded  as  a  potential  candidate  for  future  solid  electrolyte  espec¬ 
ially  in  DMFCs  [3]. 

Some  experimental  studies  have  also  been  carried  out  to  under¬ 
stand  the  microstructure  and  the  ionic  transport  inside  the  SPEEKK 
membrane.  For  example,  atomic  force  microscopy  (AFM)  and  scan¬ 
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ning  electron  microscope  (SEM)  methods  were  used  to  investigate 
the  surface  morphology  of  the  SPEEKK  membrane  [4,5].  Elec¬ 
trophoretic  nuclear  magnetic  resonance  (ENMR)  was  applied  to 
measure  the  electro-osmotic  drag  coefficients  in  Nation  and  SPEEKK 
membranes  [6].  The  states  of  ionic  aggregations  were  observed 
by  using  transmission  electron  microscope  (TEM),  small  angle 
X-ray  scattering  (SAXS),  Fourier  transform  infrared  (FTIR)  tech¬ 
niques  and  etc.  [7-9].  In  contrast  to  several  experimental  literatures 
reported  above,  only  few  theoretical  investigations  have  been  pre¬ 
sented.  Paddison  performed  ab  initio  molecular  orbital  calculations 
to  investigate  the  proton  dissociation  in  the  hydrophilic  parts  of 
hydrated  Nation  and  SPEEKK  membranes  [10];  the  proton  trans¬ 
port  phenomenon  within  the  membrane  pores  were  studied.  Using 
a  molecular  structure-function  modeling  method,  the  same  author 
computed  the  proton  friction  and  diffusion  coefficients  inside  the 
SPEEKK  membrane  at  various  levels  of  hydration  [11].  However, 
the  membrane  morphology  and  the  distribution  of  water  clusters 
in  the  SPEEKK  membrane  have  not  been  explored  in  the  previ¬ 
ous  investigations.  This  paper  intends  to  further  investigate  and 
to  compare  the  molecular  structure  and  transport  dynamics  of 
hydrated  Nafion  and  SPEEKK  membranes  by  means  of  molecular 
dynamics  technique ;  which  is  an  extension  to  our  previous  research 
[12,13].  The  morphology  difference  between  hydrated  Nafion  and 
SPEEKK  membranes  will  be  visualized  from  our  simulation  results. 
Methanol  permeability  in  these  two  membranes  will  also  be  inves¬ 
tigated.  Details  of  molecular  modeling  and  simulation  techniques 
are  described  in  the  following  section. 
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Nomenclature 

A  parameters  for  dihedral  potential  function 

(kcal  mol-1) 

a  parameters  for  12-6  potential  function  (A12) 

b  parameters  for  12-6  potential  function  (A6) 

D  diffusion  coefficient  (A2  ps-1 ) 

I<E  kinetic  energy  (kcal  mol-1 ) 

ka  force  constants  for  valence  angle  potential  function 

(kcal  mol-1  rad-2) 

kb  force  constants  for  bond  potential  function 

(kcal  mol-1  A-2) 

m  parameters  for  dihedral  potential  function 

N  number  of  atoms 

ry  distance  between  atom  i  and  atom  j  (A) 

rn  position  vectors  (A) 

r0  equilibrium  bond  length  (A) 

t  time  (ps) 

U  potential  energies  (kcal  mol-1 ) 

V  volume  (A3 ) 

8  parameters  for  dihedral  potential  function  (deg) 

0  instant  valence  angle  (rad) 

Oq  equilibrium  valence  angle  (rad) 

p  number  density  (A-3 ) 

0  instant  torsion  angle  (deg) 


s 

(a)  Nafion 


(b)  SPEEKK 

Fig.  2.  Simplified  molecular  structure  of  (a)  Nafion  and  (b)  SPEEKK  fragments  in  the 
simulation  (F:  fluorine  in  deep  blue;  C:  carbon  in  green;  H:  hydrogen  in  light  blue; 
S:  sulfur  in  yellow;  0:  oxygen  in  red).  (For  interpretation  of  the  references  to  color 
in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 


2.  Molecular  modeling  and  computer  simulation 

2 A.  Molecular  modeling 

The  molecular  simulation  system  mainly  consists  of  the  polymer 
electrolyte  (case  1  is  the  Nafion  and  case  2  is  the  SPEEKK  fragments), 
hydronium  ions,  and  solvent  molecules,  such  as  pure  water  or  water 
plus  methanol  molecules.  Fig.  1  shows  the  chemical  structure  of  the 
Nafion  and  the  SPEEKK  membrane.  For  simplification  in  molecular 
dynamics,  the  realistic  membrane  was  represented  by  fragments 
of  the  chemical  structure  instead  of  the  complete  compositions. 
The  long  backbone  of  the  Nafion  (and  the  SPEEKK)  membrane  is 
repeated  and  can  be  truncated;  the  pendant  side  chain  attached 
to  the  main  backbone  is  maintained  to  avoid  affecting  the  original 
characteristics.  Fig.  2  shows  the  simplified  Nafion  and  SPEEKK  frag- 
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(a)  Nafion  (x=5- 10,  y-1000  andz=l-2) 


ments  in  the  simulation;  the  initial  conformation  for  each  fragment 
was  constructed  by  a  semi-empirical  charge  distribution  method 
(Austin  Model  1 )  using  quantum  chemical  calculation  [14].  The  cal¬ 
culation  was  considered  to  be  converged  when  the  gradient  of  the 
atomic  force  was  less  than  the  preset  value  (0.01  kcal  mol-1  A-1  in 
our  case).  Once  convergence  was  attained,  an  optimal  geometry  of 
the  fragment  was  constructed  and  the  distribution  of  charges  was 
assigned. 

Each  proton  that  migrates  inside  the  electrolyte  is  assumed  to 
form  a  hydronium  ion  (Fl30+)  with  one  water  molecule.  The  hydro¬ 
nium  ion  is  the  simplest  type  of  an  oxonium  ion,  which  migrates 
towards  the  cathode  by  the  vehicle  and  the  Grotthuss  mechanisms 
[  15 ].  To  achieve  electrical  neutrality,  the  number  of  H3 0+  is  assigned 
to  equal  the  number  of  sulfonic-acid  (S03-)  groups.  Table  1  sum¬ 
marizes  the  simulation  conditions  of  the  investigated  system.  We 
chose  those  cases  of  water  uptakes  according  to  the  experimental 
conditions  in  Ref.  [9].  The  definition  of  water  uptake  is  the  total 
weight  of  water  divided  by  the  weight  of  the  dry  membrane  in 
experiment;  however,  in  the  simulation  we  transform  them  into  the 
ratio  of  the  number  of  molecules.  The  number  of  water  molecules 
in  the  Nafion  is  set  to  be  same  as  that  in  the  SPEEKK;  and  the  total 
number  of  water  and  methanol  molecules  is  fixed  (186  in  our  case). 
This  was  to  test  the  methanol  permeability  at  different  methanol 
concentration  in  the  later  case  in  Table  2.  Water  and  methanol 
molecules  are  initially  assigned  to  distribute  randomly  inside  the 
simulation  system. 


Ov 


SO>'H 


Fig.  1.  Chemical  structure  of  (a)  the  Nafion  and  (b)  the  SPEEKK  membrane. 


Table  1 

Summary  of  simulation  conditions. 


Membrane 


Nafion 

SPEEKK 

Water  uptake  (%) 

10.00 

9.08 

16.00 

14.53 

No.  of  fragments 

32 

32 

32 

32 

No.  of  hydronium  ions 

32 

32 

32 

32 

No.  of  water  molecules 

186 

166 

186 

166 

No.  of  methanol  molecules 

0 

20 

0 

20 

Water  uptake  (%)  = 


total  weight  of  water  molecules 
weight  of  dry  membrane 


X  100. 
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Table  2 

Density  and  volume  of  the  MD  simulation  system  obtained  from  the  1.5-ns  NPT 
simulation. 


Membrane 

Nafion 

SPEEKK 

Weight  percentage  wt%  (CH3OH) 

0 

15.11 

0 

15.11 

Molar  concentration  M  (CH3OH) 

0 

4.56 

0 

4.56 

Density  (gem-3) 

1.83 

1.76 

1.21 

1.17 

Volume  (A3) 

(34.0)3 

(34.5)3 

(59.3  )3 

(60.0)3 

wt%  (CH3OH)  = 
M  (CH3OH)  = 


total  weight  of  CH3  OH 
total  weight  ofCH3  OH,  H20  and  H30+  ' 
moleofCH3OH 


total  volume  ofCH3  OH,  H20  and  H30+  ' 


2.2.  Computer  simulation 

In  the  Hamiltonian  system,  the  total  energy  is  expressed 
by: 


^total  —  Winter  +  ^intra  +  KE  (1) 

where  Htotal  is  the  total  energy,  Dinter  is  the  inter-molecular  potential 
energy,  Dintra  is  the  intra-molecular  potential  energy,  and  I<E  is  the 
kinetic  energy.  The  inter-molecular  potential  includes  the  van  der 
Waals  potential  and  the  electrostatic  potential.  The  van  der  Waals 
force  was  expressed  by  the  typical  12-6  Lennard-Jones  function. 


For  multi-atom  molecules,  the  intra-molecular  potential  includes 
the  valence  angle  potential  (Dangle),  the  bond  potential  (Dbond),  and 
the  dihedral  angle  potential  (Ddihedrai).  These  potential  functions 
are  expressed  by: 

Uangle(0)=  y(C0S  6>-COS0o)2  (2) 

Ubond(ry)=|(ry-r0)2  (3) 

Udihedral(</>)  =  A[1  +  COS -  5)]  (4) 

where  ka,  kb  are  force  parameters  of  the  angle  and  the  bond  poten¬ 
tials;  A ,  m  and  8  are  parameters  in  the  dihedral  potential  function; 
0  and  0  are  the  bond  and  the  torsion  angles.  These  parameters 
were  evaluated  from  the  DREIDING  force  field  developed  by  Mayo 
et  al.  [16].  For  methanol  molecule,  the  three-site  potential  model 
[17]  was  employed.  The  force  parameters  of  the  O-H  and  0-CH3 
stretching  and  CH3-0-H  bending  were  obtained  from  the  all-atom 
force  fields  [18].  The  rigid  three-site  SPC/E  model  [19]  was  adopted 
for  water  molecules  and  the  force-field  parameters  of  the  hydro- 
nium  ion  were  obtained  from  the  classical  force  fields  presented  by 
Burykin  and  Warshel  [20]. 

Molecular  dynamics  (MD)  simulation  was  carried  out  on  an  IBM 
P690  workstation  using  the  software  DLPOLY  [21  ].  A  period  of  1.5-ns 
NPT  (constant  number  of  particles,  constant  pressure  and  constant 
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Fig.  3.  Total  energy  and  temperature  variation  of  the  MD  simulation  system  for  (a)  Nafion  without  methanol,  (b)  Nafion  with  methanol,  (c)  SPEEKK  without  methanol  and 
(d)  SPEEKK  with  methanol  (all  at  300 1<). 
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temperature)  ensemble  simulation  was  performed  to  tune  the  sim¬ 
ulation  system  to  reach  a  proper  density.  It  was  followed  by  another 
1.5-ns  NVT  (constant  number  of  particles,  constant  volume  and 
constant  temperature)  simulation.  All  MD  simulations  were  per¬ 
formed  at  300  K  (27  °C)  and  the  equations  of  motion  were  solved 
with  a  Verlet  scheme  [22]  with  a  time  step  1  fs.  The  Lennard -Jones 
and  electrostatic  potentials  were  truncated  at  10  A.  All  covalent 
bonds  were  kept  rigid  with  the  SHAKE  algorithm  [23]  and  three- 
dimensional  periodic  boundary  conditions  were  applied  in  the 
simulation.  The  Eward  summation  method  [24]  was  used  for  the 
calculation  of  the  electrostatic  potential  at  the  periodic  boundary 
condition. 


3.  Results  and  discussion 

3.1.  System  equilibrium 

Table  2  shows  the  calculated  density  and  volume  of  the 
equilibrated  system  obtained  from  the  preliminary  1.5-ns  NPT  sim¬ 
ulation.  This  was  to  emulate  the  membrane  swelling  with  uptake 
of  methanol.  The  weight  percentage  of  methanol  is  calculated  as 
the  ratio  between  the  total  weight  of  methanol  molecules  and  the 
total  weight  of  methanol  plus  water  molecules  plus  hydronium 
ions.  The  molar  concentration  of  methanol  is  defined  as  the  number 
of  moles  of  methanol  divided  by  the  total  volume  of  methanol  plus 
water  molecules  plus  hydronium  ions.  The  total  volume  is  evaluated 
from  the  NPT  simulation  based  on  the  input  from  Table  1.  Before 
analyzing  the  simulation  results,  we  checked  the  equilibration  of 
the  system  by  monitoring  the  total  energy  and  the  temperature 
variation  during  the  formal  1.5-ns  NVT  simulation.  Fig.  3  shows 
that  both  total  potential  energy  and  temperature  are  converged  in 
each  simulation  case.  Fluctuations  in  the  total  energy  output  are 
less  than  ±5%,  which  are  caused  by  the  exchange  of  the  potential 
and  the  kinetic  energy.  Checking  from  the  detailed  data,  the  elec¬ 
trostatic  potential  contributes  more  than  50%  of  the  total  energy, 
implying  that  the  influence  from  the  Coulombic  force  is  the  most 
significant;  while  the  influence  from  the  dihedral  potential  is  the 
least. 


3.2.  Membrane  morphology  analysis 

The  ion  conduction  properties  of  the  membrane  closely  relate 
to  the  internal  structure  such  as  the  spatial  distribution  of  ionic 
sites  in  the  membrane.  The  snapshots  of  the  membrane  mor¬ 
phology  from  the  MD  simulation  can  be  seen  in  Fig.  4(a)-Nafion 
and  (b)-SPEEKK.  Both  hydrated  Nation  and  SPEEKK  membranes 
form  the  phase  separation  of  hydrophobic  and  hydrophilic  parts. 
Our  simulation  results  show  a  good  agreement  with  the  experi¬ 
mental  observation  from  the  Ref.  [25].  The  hydrophobic  region  is 
constituted  by  the  backbone  (green  cloud)  of  the  polymer,  while 
the  hydrophilic  region  contains  water  molecules  (gray  cloud), 
hydronium  ions  (light  blue  balls)  and  sulfonic-acid  groups  (yel¬ 
low  and  red  balls).  In  comparison  with  the  hydrated  Nation  in 
Fig.  4(a),  the  hydrated  SPEEKK  in  Fig.  4(b)  forms  less  hydropho¬ 
bic/hydrophilic  phase  separations.  The  water  channels  in  the  Nation 
are  approximately  1-3  nm  in  diameter  while  those  in  the  SPEEKK 
are  approximately  0.5-1  nm.  The  nanoscopic  simulation  reveals 
that  the  water  channels  in  the  SPEEKK  are  narrower  and  less  in 
quantity  than  those  in  the  Nation.  These  features  show  a  larger 
interface  between  the  hydrophobic  and  hydrophilic  regions  in  the 
SPEEKK;  implying  that  the  average  distance  between  the  neighbor¬ 
ing  sulfonic-acid  groups  in  the  SPEEKK  is  larger  than  that  in  the 
Nation  membrane. 


(b)  SPEEKK 


Fig.  4.  Snapshots  of  the  membrane  morphology  from  the  MD  simulation  for  (a) 
Nation  and  (b)  SPEEKK  membranes  (light  blue  balls:  hydroniums;  red  balls:  oxygen 
atoms;  yellow  balls:  sulfur  atoms;  gray  cloud:  water  molecules;  green  cloud:  back¬ 
bones  of  the  polymer).  (For  interpretation  of  the  references  to  color  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 


3.3.  Molecular  structure  analysis 

The  radial  distribution  function  (RDF),  denoted  as  g(r),  is  an 
indication  of  the  local  molecular  structure  and  is  defined  mathe¬ 
matically  by 

(N(r,  Ar))  p(r)  local  number  density 
^  r  _  (l/2)NpV(r,  Ar)  -  p  ~  system  number  density 

where  ( )  indicates  a  temporal  average;  N(r,  Ar)  is  the  number  of 
atoms  within  the  spherical  shell  of  radius  between  r  and  r+  Ar;  N 
is  the  total  number  of  atoms  in  the  system;  p  is  the  number  density 
and  V(r,  Ar)  is  the  volume  of  the  shell.  Fig.  5(a)  shows  the  RDFs  of 
molecule  pairs  of  the  carbon-water  (in  which  water  is  represented 
by  the  oxygen  atom  in  it)  and  the  fluorine-water  for  the  Nation 
membrane;  there  the  RDFs  increase  almost  with  the  radius  and  the 
height  of  RDFs  is  less  than  unity,  implying  that  these  two  atoms, 
F  and  C,  are  hydrophobic  during  the  simulation  period.  The  RDF 
variation  tendency  observed  in  Fig.  5(b),  RDFs  of  molecule  pairs 
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Fig.  5.  RDFs  of  (a)  C-Oh2o,  F-Oh2o  pairs  for  Nafion  and  (b)  C-Oh2o.  H-Oh2o  pairs  for 
SPEEKK  at  300  K. 

of  carbon-water  and  hydrogen-water  for  the  SPEEKK  membrane, 
reveals  that  both  carbon  and  hydrogen  atoms  are  hydrophobic  due 
to  RDFs  stay  around  0.8  after  2.8  A  and  3.8  A,  respectively.  The  RDFs 
of  oxygen  atoms  (denoted  as  Os)  in  the  sulfonic-acid  groups  toward 
water  molecules  for  the  Nafion  and  the  SPEEKK  are  shown  in  Fig.  6. 
First  peaks  are  observed  at  the  radius  about  2.5  A  for  the  Nafion 
and  2.75  A  for  the  SPEEKK,  which  indicate  that  water  molecules 
tend  to  aggregate  near  the  sulfonic-acid  groups  in  both  membranes. 
Molecules  staying  at  the  first  peak  of  the  RDF  diagram  correspond 
to  the  “bound  state”  and  are  restrained  from  moving  freely.  In  addi¬ 
tion,  the  height  of  the  first  peak  of  the  Nafion  is  greater  than  that 
of  the  SPEEKK,  implying  that  the  local  number  density  of  water 


Fig.  6.  RDFs  of  Os-Oh2o  pair  for  Nafion  and  SPEEKK  membranes  at  300 1<. 


Fig.  7.  RDFs  of  (a)  Oh2o-Oh2o  and  (b)  Oh2o-0H3o+  pairs  for  Nafion  and  SPEEKK 
membranes  at  300 1<. 

molecules  associated  with  the  sulfonic-acid  groups  in  the  Nafion  is 
greater  than  that  in  the  SPEEKK.  Fig.  7(a)  and  (b)  shows  the  RDFs  of 
molecule  pairs  of  water-water  and  water-hydronium  for  the  Nafion 
and  the  SPEEKK  membranes.  Significant  first  peaks  in  both  (a)  and 
(b)  indicate  that  water  molecules  tend  to  form  clusters  and  they 
can  drag  the  hydronium  ions  to  travel  in  these  two  membranes.  The 
dragging  effect  is  less  significant  in  the  SPEEKK  than  in  the  Nafion. 

Fig.  8  shows  the  distributions  of  water  clusters  and  their  sizes 
in  the  Nafion  and  the  SPEEKK  membranes.  Water  molecules  are 
assumed  to  belong  to  the  same  cluster  if  the  distance  between  the 
oxygen  atoms  from  different  water  molecules  is  less  than  3.5  A.  It 
is  shown  that  water  molecules  in  the  SPEEKK  form  more  clusters 
in  comparison  with  the  Nafion,  54  versus  58  as  indicated  by  the 
arrows  in  the  diagram;  while  the  average  size  of  water  clusters  in 
the  SPEEKK  is  smaller  than  that  in  the  Nafion.  Water  molecules 
are  considered  to  form  a  hydrophilic  aggregation  with  the  sulfonic- 
acid  groups  when  the  distance  between  the  oxygen  atom  of  the 
RS03_  and  the  oxygen  atom  of  the  water  is  within  a  preset  value 
(2.5  A  for  the  Nafion  and  2.75  A  for  the  SPEEKK,  respectively).  The 
chosen  value  corresponds  to  the  average  distance  of  the  first  hydra¬ 
tion  shell  around  the  sulfonic-acid  groups.  The  average  numbers  of 
water  molecules  around  a  specific  sulfonic-acid  group  are  4.35  for 
the  Nafion  and  3.07  for  the  SPEEKK.  The  results  of  water  clusters 
distributions  and  hydrophilic  aggregations  analysis  agree  well  with 
the  membrane  morphology  shown  in  Fig.  4,  where  the  hydrated 
Nafion  forms  a  larger  hydrophobic/hydrophilic  phase  separation 
and  the  water  channels  in  the  Nafion  are  wider  than  those  in  the 
SPEEKK.  These  features  lead  to  a  rise  in  the  size  of  water  clusters 
and  the  hydrophilic  aggregation  in  the  Nafion  than  in  the  SPEEKK 
membrane. 
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Fig.  8.  Distribution  of  water  molecule  number  versus  water  cluster  number  inside 
(a)  Nation  and  (b)  SPEEKK  membranes. 


3.4.  Molecular  mobility  analysis 

The  mobility  of  molecules  can  be  analyzed  using  the  mean 
square  displacement  (MSD),  which  is  defined  mathematically  by: 

N 

MSD^E([^(to+t)-^(to)l2)  (6) 

n=l 

where  N  is  the  number  of  molecules;  rn(t)  is  the  position  vector 
of  molecule  n  at  time  t  and  to  is  the  initial  time  step.  Fig.  9  shows 
the  MSDs  of  hydronium  ions  and  water  molecules  for  the  Nation 
and  the  SPEEKK  membranes.  Linear  tendencies  of  the  MSD  curves 
infer  that  the  hydronium  ions  and  water  molecules  are  diffusing 
continuously  in  the  system  during  the  simulation.  Both  hydronium 
ions  and  water  molecules  show  a  greater  slope  of  MSD  curves  in  the 
Nation  than  in  the  SPEEKK.  Diffusion  coefficients  can  be  evaluated 
from  the  Einstein  relation  as  below 

D=^t1™lE^(to+t)-^(to)2>  (7) 

n= 1 

where  D  is  the  diffusion  coefficient,  N  is  the  number  of  molecules 
or  ions,  and  rn(t)  is  the  position  vector  of  molecule  or  ion  n  at  time 
t.  Table  3  shows  the  comparison  of  diffusion  coefficients  of  water 
molecules  in  the  Nation  and  the  SPEEKK  between  simulations  and 
experiments.  Our  predictions  are  in  reasonable  agreement  with  the 
experimental  results  from  Refs.  [9]  and  [12]  (0.42  x  10-5  cm2  s-1  for 
the  Nation  and  0.30  x  10-5  cm2  s-1  for  the  SPEEKK,  respectively). 
MD  predictions  of  the  ionic  conductivity  of  hydronium  ions  in 
the  Nation  and  in  the  SPEEKK  membranes  are  0.011  S  cm-1  and 


(a)  500  r 


500  1000 

Time  (ps) 


1500 


Fig.  9.  MSDs  of  (a)  water  molecules  and  (b)  hydronium  ions  for  Nation  and  SPEEKK 
membranes  at  300 1<. 


Table  3 

Comparison  of  diffusion  coefficient  of  water  molecules  in  the  Nation  and  the  SPEEKK 
between  simulation  results  and  similar  experiments  [9,12]. 


Diffusion  coefficient  (cm2  s_1 ) 

Membrane 

of  water  molecules 

Nation 

SPEEKK 

MD  prediction 

0.55  x  10“5 

0.39  x  10“5 

Experiment 

0.42  x  10-5 

0.30  x  10-5 

0.003  S  cm-1,  respectively.  They  are  about  one  third  of  the  exper¬ 
imental  results  reported  from  Refs.  [8,26]  as  shown  in  Table  4. 
The  discrepancy  comes  from  the  molecular  model  that  we  only 
consider  the  vehicle  mechanism  for  modeling  the  hydronium  trans¬ 
port;  while  the  transfer  of  proton  hopping  (i.e.,  the  Grotthuss 
mechanism)  between  hydronium  ions  and  water  molecules  is  not 
included.  The  later  mechanism  needs  more  detailed  quantum  sim¬ 
ulation.  Summarily,  the  mobilities  of  hydronium  ions  and  water 
molecules  in  the  SPEEKK  are  lower  than  those  in  the  Nation.  This  can 
be  explained  by  the  narrower  water  channels  and  less  hydrophilic 
aggregations  in  the  SPEEKK  restrict  the  migration  of  hydronium 
ions  and  water  molecules. 


Table  4 

Comparison  of  ionic  conductivity  of  hydronium  ions  in  the  Nation  and  the  SPEEKK 
between  simulation  results  and  similar  experiments  [8,26]. 


Ionic  conductivity  (S  cm-1 )  of  hydronium  ions 

Membrane 

Nation 

SPEEKK 

MD  prediction 

0.011 

0.003 

Experiment 

0.038 

0.013 
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Fig.  10.  MSDs  of  methanol  permeability  in  Nation  and  SPEEKK  membranes  at  300 1<. 


Table  5 

Comparison  of  diffusion  coefficient  of  methanol  molecules  in  the  Nafion  and  the 
SPEEKK  between  simulation  results  and  similar  experiments  [8,27]. 


Diffusion  coefficient  (cm2  s-1 ) 

Membrane 

of  methanol  molecules 

Nafion 

SPEEKK 

MD  prediction 

0.55  x  10-5 

0.024  x  10“5 

Experiment 

0.50  x  10“5 

0.017  x  10“5 

3.5.  Methanol  permeability  analysis 

Methanol  permeation  through  the  membrane  deteriorates  the 
performance  of  DMFCs  significantly.  This  can  be  analyzed  from 
the  MD  simulation  results.  Fig.  10  shows  the  MSDs  of  methanol 
molecules  for  the  Nafion  and  the  SPEEKK  membranes.  The  slope 
of  the  MSD  trace  in  the  Nafion  shows  much  steeper  than  in  the 
SPEEKK;  implying  that  the  methanol  permeation  will  be  more  seri¬ 
ous  in  the  Nafion  membrane.  Table  5  shows  the  comparison  of 
methanol  diffusion  coefficients  in  the  Nafion  and  in  the  SPEEKK 
between  our  simulation  and  experiments  [8,27].  They  are  in  rea¬ 
sonable  agreement.  The  diffusion  coefficient  of  methanol  in  the 
SPEEKK  is  much  lower  than  that  in  the  Nafion.  A  rationalization 
of  this  result  is  attributed  to  (a)  the  sulfonic-acid  groups  are  dis¬ 
persed  throughout  the  SPEEKK  membrane;  (b)the  phase  separation 
in  the  SPEEKK  is  less  significant;  and  (c)  the  water  transport  chan¬ 
nels  in  SPEEKK  are  narrower  compared  with  those  in  the  Nafion.  We 
can  draw  a  conclusion  here  that  methanol  molecules  travel  with 
water  molecules  and  hydronium  ions  via  the  hydrophilic  sulfonic- 
acid  groups.  Narrower  water  channels  in  the  SPEEKK  make  the 
methanol  diffusion  coefficient  in  the  SPEEKK  much  lower  than  that 
in  the  Nafion.  The  results  indicate  that  the  SPEEKK  membrane  is 
able  to  reduce  the  methanol  permeability  from  the  Nafion  mem¬ 
brane  without  much  sacrifice  of  water  and  hydronium  diffusion  in 
DMFCs. 

4.  Conclusions 

The  molecular  structure  and  transport  dynamics  of  hydrated 
Nafion  and  SPEEKK  membranes  are  investigated  in  this  paper  by 


means  of  molecular  simulation  techniques.  It  is  found  that  both 
hydrated  Nafion  and  SPEEKK  membranes  form  the  phase  separa¬ 
tion  of  hydrophobic  and  hydrophilic  regions;  the  hydrated  SPEEKK 
forms  less  phase  separations  in  comparison  with  the  Nafion.  Nar¬ 
rower  and  more  branched  water  channels  in  the  SPEEKK  result  in  a 
larger  interface  between  the  hydrophobic  and  hydrophilic  regions. 
These  characteristics  lead  to  a  rise  in  the  number  of  water  clusters 
and  a  decrease  of  water  molecule  number  in  each  cluster  in  the 
SPEEKK  membrane.  Narrower  water  channels  in  the  SPEEKK  make 
the  methanol  diffusion  coefficient  in  the  SPEEKK  much  lower  than 
that  in  the  Nafion.  The  results  indicate  that  the  SPEEKK  membrane 
is  able  to  reduce  the  methanol  permeability  from  the  Nafion  mem¬ 
brane  without  much  sacrifice  of  water  and  hydronium  diffusion  in 
DMFCs.  This  paper  has  established  a  molecular  simulation  tech¬ 
nique  to  investigate  the  nano-technology  to  further  improve  the 
methanol  permeation  problem  in  DMFCs  by  means  of  designing 
new  molecular  structure  in  proton-conducting  electrolytic  mem¬ 
branes.  This  simulation  technique  is  an  effective  and  efficient  way 
to  help  designing  new  materials  in  the  future. 
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